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Abstract 

Gravitino dark matter, together with thermal leptogenesis, implies an upper bound 
on the masses of superparticles. In the case of broken R-parity the constraints from 
primordial nucleosynthesis are naturally satisfied and decaying gravitinos lead to char- 
acteristic signatures in high energy cosmic rays. We analyse the implications for su- 
pergravity models with universal boundary conditions at the grand unification scale. 
Together with low-energy observables one obtains a window of superparticle masses, 
which will soon be probed at the LHC, and a range of allowed reheating temperatures. 



1 Introduction 



Standard thermal leptogenesis [1] provides a simple and elegant explanation of the origin of 
matter. It is a natural consequence of the seesaw mechanism, and it is perfectly consistent 
with the small neutrino masses inferred from neutrino oscillation data [2]. 

Thermal leptogenesis works without and with supersymmetry. In the latter case, how- 
ever, there is a clash with the 'gravitino problem' [3-5]: the large temperature required 
by leptogenesis exceeds the upper bound on the reheating temperature from primordial 
nucleosynthesis (BBN) in typical supergravity models with a neutralino as lightest super- 
particle (LSP) and an unstable gravitino. If the gravitino is the LSP, the condition that 
relic gravitinos do not overdose the universe yields an upper bound on the reheating tem- 
perature [6]. Furthermore, the next-to-lightest superparticle (NLSP) is long lived, and one 
has to worry about the effect of NLSP decays on nucleosynthesis. 

It is remarkable that, despite these potential problems, a large leptogenesis tempera- 
ture of order 10 10 GeV can account for the observed cold dark matter in terms of thermally 
produced relic gravitinos [7]. Requiring consistency with nucleosynthesis yields constraints 
on the superparticle mass spectrum. Due to improved analyses of BBN, the original pro- 
posal of a higgsino NLSP is no longer viable, and also other possible NLSPs are strongly 
constrained. The case of a stau NLSP is cornered by bounds following from catalyzed 
production of 6 Li [8], with the possible exception of a large left-right mixing in the stau 
sector [9]. In some models a sneutrino [10] or a stop [11] can still be a viable NLSP. 

Recently, it has been shown that in the case of small R-parity and lepton number break- 
ing, such that the baryon asymmetry is not erased by sphaleron processes [12], thermal 
leptogenesis, gravitino dark matter and primordial nucleosynthesis are naturally consis- 
tent [13]. Although the gravitino is no longer stable, its decay into standard model (SM) 
particles is doubly suppressed by the Planck mass and the small R-parity breaking pa- 
rameter. Hence, its lifetime exceeds the age of the universe by many orders of magnitude, 
and it remains a viable dark matter candidate [14]. Gravitino decays lead to characteris- 
tic signatures in high energy cosmic rays. The produced flux of gamma-rays [13-17] and 
positrons [17,18] may explain the observed excess in the EGRET [19] and HEAT [20] data. 
This hypothesis will soon be tested by the satellite experiments FGST and PAMELA. 

In this paper we study the implications of leptogenesis and gravitino dark matter with 
broken R-parity on the mass spectrum of superparticles. Since the unification of gauge 
couplings in the minimal supersymmetric extension of the standard model (MSSM) is one 
of the main motivations for low-energy supersymmetry, we shall focus on versions of the 
MSSM with universal boundary conditions for scalar and gaugino masses at the grand 
unification (GUT) scale. As we shall see, the corresponding spectrum of superparticle 
masses will be fully covered at the LHC. This is the main result of our analysis. 

After some comments on R-parity violation in Section 2, we discuss the lower bound on 
the reheating temperature from leptogenesis and the upper bound on the NLSP mass from 
gravitino dark matter in Section 3. Section 4 deals with constraints on MSSM parameters 
from low-energy observables, and the results of our numerical analysis are presented in 
Section 5, followed by some conlusions in Section 6. 
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2 Constraints on R-parity violation 



Phenomenological aspects of R-parity violation have been widely discussed in the literature 
[21]. Here we are interested in the case of small R-parity and lepton number breaking which 
was investigated in [13, 14, 17]. The details strongly depend on the flavour structure of R- 
parity violating couplings and the pattern of supersymmetry breaking. For completeness, 
we recall in the following the order of magnitude of bounds on R-parity violating couplings, 
the corresponding lifetimes of gravitino and NLSP, and in particular the dependence on 
the gravitino mass. 

Stringent constraints on the lepton number and R-parity violating interactions 

W A l=i = \ikjhe c jl k + KjidiQjh (2.1) 

are imposed by baryogenesis. Both operators contain lepton doublets. Together with 
sphaleron processes they therefore influence the baryon asymmetry at high temperature 
in the early universe. The requirement that an existing baryon asymmetry is not erased 
before the electroweak transition typically implies [12] 

A,A'<10~ 7 . (2.2) 

Remarkably, for such a small breaking of R-parity a gravitino LSP has a lifetime much 
longer than the age of the universe [14] because of the double suppression of the decay rate 
by the inverse Planck mass and the R-parity breaking coupling. One then obtains for the 
gravitino lifetime (cf. [13]) 

Ts/2 , w « s ^y\^ { ^ r \ ( , 3) 

where in ~ (9(100 GeV) is a characteristic supersymmetry breaking mass scale. In the 
case of light gravitinos, m 3 / 2 <C m, where only the decay into photon neutrino pairs is 
kinematically allowed, i] = 1 has been assumed in [13]. For heavier gravitinos, decays 
into W-boson lepton and Z-boson lepton pairs are also possible, and we only know that 
7] — 0(1) [17]. In particular, the relation between gravitino lifetime and gravitino mass 
depends on the pattern of supersymmetry breaking. 

In the case of a small breaking of R-parity, with an unstable gravitino LSP, the NLSP 
lifetime becomes very short, 

For couplings A, A' > 10~ 14 , the NLSP lifetime becomes shorter than 10 3 s. In case of 
a stau NLSP, superparticle decays then do not affect the primordial abundances of light 
elements. Hence, baryogenesis, primordial nucleosynthesis and gravitino dark matter can 
be consistent in the range 

10- 14 < A, A' < 10~ 7 . (2.5) 
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For a bino NLSP, a lifetime shorter than 0.1 s, i.e., couplings A, A' > 10~ 12 are required by 
consistency with BBN. 

The analysis of constraints on the superpotential terms (12.11) can be extended to general 
R-parity breaking mass terms [17], yielding again a range of allowed parameters. One 
finds that possible contributions to neutrino masses are negligable, once the cosmological 
constraints are satisfied. 

Decaying gravitino dark matter can contribute to the EGRET and HEAT anomalies 
for a gravitino lifetime r 3 / 2 ~ 10 26 s. For a gravitino mass m 3 / 2 ~ 10 GeV, and assuming 
i] ~ 1 in Eq. (12.31) . this requires R-parity violating couplings A ~ 10~ 7 . As we shall see, 
universal boundary conditions for gaugino masses favour larger gravitino masses, in the 
range 

10 GeV < m 3/2 < 500 GeV , (2.6) 

which, for fixed gravitino lifetime and 77 ~ 1, corresponds to the range of R-parity violating 
couplings 

10~ 10 < A < 10" 7 . (2.7) 

Note that for couplings below ~ 10~ 9 , most NLSPs decay outside the detector. However, 
for couplings above ~ 10~ n , corresponding to lifetimes shorter than ~ 10~ 3 s, some NLSP 
decays may still be observable in the detector [22]. 

How can the phenomenologically required small R-parity violating couplings arise? 
In [13] an example was presented, where the spontaneous breaking of R-parity is tied to 
B-L breaking. Recently, it has been shown that also the breaking of left-right symmetry 
can lead to small R-parity breaking [23]. 



3 Thermal leptogenesis 



Let us now consider standard thermal leptogenesis as the source of the cosmological baryon 
asymmetry. In the high-temperature phase of the early universe thermally produced right- 
handed neutrinos generate an asymmetry in B-L, which leads to a baryon asymmetry via 
sphaleron processes. In the case of hierarchical right-handed neutrinos, and neglecting 
flavour effects, the baryon density relative to the photon density is given by (cf. [2]) 



— ~ -1.04 x 10~ 2 eiK, 

n 7 



(3.1) 



where e\ is the CP asymmetry in the decay of the lightest right-handed neutrino N\ into a 
pair of lepton (L) and Higgs (H u ) doublets, the efficiency factor k represents the effects of 
washout and scattering processes, and we have assumed a supersymmetric thermal plasma. 
The CP asymmetry ei satisfies an upper bound because of the seesaw relation, which for 
supersymmetric models reads [24-26], 



e i 



r(AT 1 _> L + H u ) + T(N X -> L c + HI 



< 



3Mi Am 



at in 



8tt(H u } 2 mi + m 3 ' 



(3.2) 
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Here mi, with mi < m 2 < m 3 , are the mass eigenvalues of the light neutrinos and Mi is the 
mass of the right-handed neutrino Ni. The atmospheric neutrino mass squared difference is 
determined from neutrino oscillation experiments as Am^ tm ~ (2.5 ±0.2) x 10 _3 eV 2 . Note 
that the upper bound on |ei|, and therefore the maximally generated baryon asymmetry, 
increases proportional to the heavy Majorana mass M%. 

The efficiency factor k has to be determined by solving the Boltzmann equations. In 
the most interesting case of zero initial abundance of the right-handed neutrinos one finds 
for its maximal value, with and without supersymmetry, k ~ 0.2 [27,28]. Using (13 .21) . one 
then obtains from the observed baryon asymmetry [29], 

— = (6.21 ±0.16) x 10- 10 , (3.3) 

the lower bound on the right-handed neutrino mass 

Ml > 1.4 x ID" GeV (±§±f (3.4) 

at the 3d level of and Am^ tm . The corresponding lower bound on the reheating 

temperature is about a factor two smaller [30]. In the following analysis we shall therefore 
use as an estimate 

Tr > 1 x 10 9 GeV. (3.5) 

Note that this bound on the reheating temperature only applies for hierarchical right- 
handed neutrinos. In the case of quasi-degenerate heavy neutrinos it is relaxed. The 
bound also assumes thermal equilibrium, and it is modified once the reheating process is 
taking into account. For instance, in the case of reheating by inflaton decays, the bound 
increases by about a factor of two [28]. 

Relic gravitinos with masses larger than 1 GeV contribute to cold dark matter. In 
the following analysis we identify the thermally produced abundance £^3/2 ^ 2 with the 2a 
upper bound on the dark matter abundance deduced from the CMB anisotropies. From 
the WMAP 5-year results one obtains [29], 

tt 3/2 h 2 = tt DM h 2 ~ 0.1223 . (3.6) 

The thermal production of gravitinos is dominated by QCD processes. To leading order 
in the gauge coupling we find 

where m g i uino is the physical gluino mass. Note that the coefficient!]] is about a factor two 
larger than in the analysis [31]. This is due to the 2-loop running of the gluino mass, which 
has been taken into account. Electroweak contributions to thermal gravitino production 



1 Varying superparticle masses (cf. Section 4), the value can change by about 10%. 



5 



further increase the abundance by about 20%. In our numerical analysis we shall take 
this into account following [32]. Note that the gravitino production rate has an 0(1) 
uncertainty due to unknown higher order contributions and nonperturbative effects [31]. 
Resummation of thermal masses increases the production rate by about a factor of two 
[33]. We also neglect nonthermal contributions to gravitino production, in particular from 
inflaton decay [34], which are usually subdominant at the considered high temperatures. 

Our main interest are constraints on gluino and NLSP masses for gravitino dark matter. 
It is then convenient to rewrite (13.7ft as 



where the ratio £ is fixed by the boundary conditions of the soft supersymmetry breaking 
parameters. For each gravitino mass and reheating temperature, Eq. (I3.8P then gives the 
NLSP mass for which the observed dark matter density is obtained. The maximal NLSP 
mass is reached for m 3 / 2 = m-NLSP, 



In this paper, we focus on thermally produced gravitino dark matter. A high reheating 
temperature can also be consistent with leptogenesis in the case of very heavy gravitinos, 
as in anomaly mediation [35] or mirage mediation [36,37]. In those models, the gravitino 
can have a mass of about 100 TeV and thus decays before BBN starts. However, these 
models have several intrinsic difficulties. In the case of anomaly mediation, it is difficult to 
explain the g — 2 anomaly together with the b — > 57 constraint, since the gaugino masses 
are controlled by the beta functions. In mirage mediation models, one often has a light 
modulus field whose decay produces too many gravitinos [38]. Hence, the heavy gravitino 
scenario appears to be phenomenologically disfavoured. 

4 Models and low-energy observables 

In order to illustrate the implications of leptogenesis and gravitino dark matter on su- 
perparticle masses, we now study two typical boundary conditions for the supersymmetry 
breaking parameters of the MSSM at the grand unification (GUT) scale: 

(A) mo — mi/2, a = 0, tan/3, (4.1) 

with equal universal scalar and gaugino masses, mo and mi/2, respectively; in this case a 
bino-like neutralino becomes the NLSP. The second boundary condition is 

(B) m = 0, mi/2, ao = 0, tan/?, (4.2) 

which yields the right-handed stau as NLSP. In both cases, the trilinear scalar coupling ao 
is put to zero for simplicity. The ratio tan j3 of the Higgs vacuum expectation values and the 




(3.8) 




(3.9) 
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universal gaugino mass m\/2 are the two remaining independent variables. Superparticle 
masses at the electroweak scale are obtained by solving the renormalization group equations 
at 2-loop accuracy by means of SOFTSUSY 2.0.18 [39]. 

Low-energy observables yield a lower bound on superparticle masses. Since the thermal 
gravitino abundance (13 .7p increases quadratically with the gluino mass, this implies an up- 
per bound on the reheating temperature. Together with the lower bound from leptogenesis 
one then obtains a range of allowed reheating temperatures. In the same way, leptogenesis 
and gravitino dark matter yield an upper bound on superparticle masses. Combined with 
low-energy constraints, a window of allowed superparticle masses is obtained. 

One of the strongest constraints on the MSSM parameter space follows from the lower 
bound on the Higgs boson mass by LEP [40], 

m h > 114.4 GeV (95%C.L.) . (4.3) 

The bound is satisfied by enhancing radiative corrections to the Higgs potential, which 
requires a large stop mass. The parameters of the stop sector are essentially controlled by 
the gluino mass, i.e. mi/2, via the renormalization group evolutions; they are less sensitive 
to the scalar mass m . The potential is also affected by the trilinear stop coupling A t for 
sufficiently large a$. Although we put ao = in the numerical analysis, we shall comment 
on the case a ^ 0. In our analysis we use the top quark mass m t = 172.6 GeV [40]. 
Radiative corrections are taken into account at the 2-loop level by means of FeynHiggs 
2.6.4 [41]. 

When the superparticles are light, they contribute significantly to rare processes. The 
measured branching ratio Br(£> rf — > Xg'j) agrees with the SM prediction. The SUSY 
contributions are dominated by the top-charged Higgs and stop-chargino diagrams. The 
latter is enhanced by large tan j3 and interferes with the former. In our analysis we choose 
the sign of the supersymmetric Higgs mass parameter hh such that the effect of the SUSY 
contributions is reduced. Taking into account the theoretical uncertainties, we require for 
the full MSSM prediction the conservative upper and lower bounds, 

2 x 10~ 4 < Bi{B d X sl ) < 4 x 10 -4 . (4.4) 

The numerical analysis is based on SusyBSG 1.1.2 which takes NNLO contributions partially 
into account [42]. 

The two observables discussed above constrain the MSSM parameters. In contrast, the 
apparent discrepancy between the measured value of the muon anomalous magnetic mo- 
ment [43] and the SM prediction may be an effect of supersymmetry, which then favours 
a certain range of MSSM parameters. Recently, the hadronic contribution to the SM pre- 
diction has been updated using e + e~ data [44]. The current discrepancy with experiment 
is given by [45] 

a M (exp) - a^(SM) = 302(88) x 10~ n , (4.5) 

which corresponds to a 3.4<r deviation. An explanation of this discrepancy by hypothetical 
errors in the determination of the hadronic SM contribution appears unlikely [45]. In con- 
trast, supersymmetry can easily account for the discrepancy [46]. The SUSY contribution 
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is proportional to tan/3 and depends on sgn(//#). It is remarkable that the deviation from 
the SM prediciton for a M and the agreement for Br(B d — > X^) require the same sgn(yU^) in 
the case of universal gaugino masses at the GUT scale. In the following, we use FeynHiggs 
to evaluate the SUSY contribution to a M at the 2-loop level. 

Finally, the absense of pair production of heavy charged particles at LEP implies the 
approximate lower mass bound [40] 

m charged > 100 GeV . (4.6) 

In the next section we shall use superparticle masses obtained by means of SOFTSUSY. 



5 Numerical analysis 

We are now ready to determine the superparticle mass window and the allowed range of 
reheating temperatures for the two examples of universal boundary conditions at the GUT 
scale, which were discussed in the previous section. 

In Fig.CQthe upper bound (13.90 on the NLSP masses is shown for reheating temperatures 
Tr > 1 x 10 9 GeV, which is the lower bound required by leptogenesis. In case (A) with 
bino NLSP, the ratio £ = ^NLSp/^giuino, and therefore the upper bound on m NLSP , are 
essentially independent of tan j3. In contrast, for (B) with stau NLSP, one has a strong 
dependence on tan f3. The lower bound on m NL sp is determined by Br(£> d — > X s j) and the 
Higgs mass bound in case (A), and the charged particle and Higgs mass bounds in case 
(B), respectively. We find the allowed mass ranges 

(A) 130 GeV < m bino < 620 GeV , (B) 100 GeV < m stau < 490 GeV . (5.1) 

Note that in case (B) upper and lower bounds correspond to different values of tan f3. The 
muon g-2 anomaly favours small NLSP masses in the range from 100 GeV to 300 GeV. 
One also obtains upper bounds on the gravitino mass, 

(A) m 3/2 < 620 GeV , (B) m 3/2 < 490 GeV . (5.2) 

Both boundary conditions have a = 0. For negative ao, the Higgs boson potential is 
modified in such a way that the dashed line in Fig. [Q moves to the left. We have checked 
that the reheating temperature can then reach 6 x 10 9 GeV, whereas other observables are 
not much affected. We therefore obtain for the range of reheating temperatures consistent 
with leptogenesis and gravitino dark matter 

Tjt = (1 — 6) x 10 9 GeV . (5.3) 

Note that according to FeynHiggs, the theoretical uncertainty of the Higgs boson is about 
1 GeV for m h ~ 115 GeV. This corresponds to an uncertainty of 10 — 20% for the upper 
bound on the reheating temperature. 
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Figure 1: Contours of constant reheating temperature, T R — (1 — 4) X 10 9 GeV, with 
^3/2 = ^dm (solid lines) (cf. Eq. (13.61) ). The panels (a) and (b) correspond to the GUT 
boundary conditions (A) and (B) with bino-like NLSP and stau NLSP, respectively. The 
choice m 3 / 2 = 777 NL sp maximizes the reheating temperature. The gray region is excluded 
by constraints from low-energy experiments: the lower tan (3 part (left of the dashed line) 
does not satisfy the LEP Higgs mass bound; the higher tan (3 part in (a) (left of the dotted 
line) is ruled out by Br(5 (i — > X s j); the higher tan (3 part in (b) (left of the dot-dashed line) 
does not satisfy the lower mass bound on charged particles from LEP. Thermal leptogenesis 
is possible in the yellow and orange regions; the orange region is favored by the muon g — 2 
anomaly at the 2a level. 



We can also study superparticle masses as function of gravitino mass and reheating 
temperature using Eq. (13.81) . The allowed NLSP mass range then depends on tan/3. In the 
case of bino NLSP, consider as an example 

(A) tan/3 = 30, £ = = 0.17 - 0.19 . (5.4) 

^gluino 

The left panel of Fig. [2] shows the bino mass yielding the observed dark matter abundance 
as function of the gravitino mass for different reheating temperatures; the right panel is the 
corresponding plot for the gluino mass. Upper mass bounds are obtained for the smallest 
temperature of 1 x 10 9 GeV and the largest gravitino mass 777.3/2 = raunoi 

(A) tan (3 = 30 : m bino < 620 GeV , 777 gluino < 3.1 TeV . (5.5) 

For smaller gravitino masses the bounds become more stringent. For instance, for 7773/2 = 
100 GeV, one obtains 

^bino & 270 GeV , 777 g i uino < 1.5 TeV . (5.6) 
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Figure 2: Contours of constant reheating temperature in the mbi no — m 3/2 plane (a) and 
the m g i u ; no — m 3 / 2 plane (b) for boundary condition (A) with bino NLSP (see caption of 
Fig. CD for details). In the dark gray region, the gravitino is not the LSP. 



Note that these bounds are essentially independent of m and tan/3, als long as m ~ m 3 / 2 . 

In the case of stau NLSP, there is a strong dependence on tan/3. As an example, we 
consider 

(B) tan/3 =10, g= TOstau =0.16-0.17. (5.7) 

Mgluino 

£ decreases with increasing tan/3. Stau and gluino masses are shown in Fig. 3. Since the 
ratio of NLSP and gluino mass is smaller, the mass bounds are now more stringent, 

(B) tan/3 = 10: m stau < 490 GeV , m gluino < 2.8 TeV . (5.8) 

For a gravitino mass m 3 / 2 = 100 GeV, one obtains 

m stau < 240 GeV , m gluino < 1.5 TeV . (5.9) 

Let us emphasize again the effect of the theoretical uncertainty in the evaluation of 
the gravitino abundance, which is expected to be 0(1) [31]. For instance, if the gravitino 
production rate is larger by a factor 2, as suggested in [33], all reheating temperatures in 
Figs. [Tj, [2] and [3] are by a factor 2 smaller. Hence, the superparticle mass range consis- 
tent with thermal leptogenesis becomes narrower. On the other hand, a smaller gravitino 
production rate would enlarge the parameter range consistent with leptogenesis. 

Finally, let us comment on other boundary conditions. We have chosen universal gaug- 
ino masses, with m = myz or m = at the GUT scale. However, even for non-universal 
gaugino masses we obtain almost the same results. The reason is that all the bounds are 
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Figure 3: Contours of constant reheating temperature in the ?Ti stau — m 3 / 2 plane (a) and 
the m g i uin0 — m 3 / 2 plane (b) for boundary condition (B) with stau NLSP (see caption of 
Fig. Q] for details). In the dark gray region, the gravitino is not the LSP. 

controlled by the gluino mass. Reducing the gluino mass, the dark matter bound on the 
reheating temperature is relaxed, but the low-energy constraints become severer: super- 
symmetric contributions to the Higgs boson mass are suppressed, while they are enhanced 
for Br(Bd — > X s j). As a consequence, the maximal reheating temperature remains almost 
the same as in the case of universal gaugino masses. On the other hand, the low-energy 
constraints become weaker for scalar masses much larger than m^. One can then reach 
reheating temperatures ~ 10 10 GeV. 

6 Conclusions and outlook 

We have studied the implications of thermal leptogenesis and gravitino dark matter for 
the mass spectrum of superparticles. In the case of broken R-parity the constraints from 
nucleosynthesis are naturally fulfilled, and universal gaugino masses at the GUT scale are 
possible, contrary to the case of stable gravitinos. 

As an illustration, we have considered two boundary conditions which lead to a bino-like 
NLSP and a stau NLSP, respectively. Low-energy observables and gravitino dark matter 
together with thermal leptogenesis yield upper and lower bounds on NLSP and gluino 
masses, which in both cases lie within the discovery range of the LHC. It is encouraging 
that the supersymmetric explanation of the muon g — 2 anomaly favours smaller masses 
within these mass windows. 

A cosmology with leptogenesis and gravitino dark matter also leads to the prediction 
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of a maximal temperature in the early universe. In the case of universal gaugino masses 
at the unification scale we find the upper bound T™ x ~ 6 x 10 9 GeV, which is somewhat 
relaxed for large scalar masses. This bound has been obtained under the assumption of 
thermal equilibrium, which appears unlikely for a maximal temperature. Nevertheless, it is 
intriguing that the temperature T^ ax is of the same order of magnitude as the critical for the 
destabilization of compact dimensions in higher-dimensional supersymmetric theories [47]. 
The effect of the reheating process on the stabilization of extra dimensions and the relation 
to baryogenesis and dark matter require futher investigations. 

Gravitino decays produce a flux of photons and positrons, which can significantly con- 
tribute to the EGRET and HEAT anomalies for a lifetime r 3 / 2 ~ 10 26 s. If these anomalies 
are indeed related to gravitino decays, the satellite experiments FGST and PAMELA 
should soon detect characteristic features in the photon and positron spectrum, respec- 
tively. Observation of a line in the gamma-ray spectrum by FGST and a rise with sharp 
cutoff in the positron spectrum by PAMELA would lead to a determination of the grav- 
itino mass. This would considerably tighten the predictions for superparticle mass windows 
which will be probed at the LHC. 

References 

[1] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986). 

[2] For a review and references, see W. Buchmuller, R. D. Peccei and T. Yanagida, Ann. 
Rev. Nucl. Part. Sci. 55 (2005) 311; S. Davidson, E. Nardi and Y. Nir. larXiv:0802.2962l 
[hep-ph]. 

[3] S. Weinberg, Phys. Rev. Lett. 48 (1982) 1303. 

[4] J. R. Ellis, D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B 259 (1985) 175. 

[5] M. Kawasaki, K. Kohri and T. Moroi, Phys. Lett. B 625, 7 (2005); Phys. Rev. D 71, 
083502 (2005); K. Jedamzik, Phys. Rev. D 74, 103509 (2006). 

[6] T. Moroi, H. Murayama and M. Yamaguchi, Phys. Lett. B 303, 289 (1993). 

[7] M. Bolz, W. Buchmuller and M. Plumacher, Phys. Lett. B 443 (1998) 209. 

[8] M. Pospelov, Phys. Rev. Lett. 98, 231301 (2007). 

[9] M. Ratz, K. Schmidt-Hoberg and M. W. Winkler. larXiv:0808.0 829 [hep-ph]; J. Pradler 
and F. D. Steffen, arXiv:0808.2462l [hep-ph]. 

[10] T. Kanzaki, M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 75, 025011 (2007). 

[11] J. L. Diaz-Cruz, J. R. Ellis, K. A. Olive and Y. Santoso, JHEP 0705 (2007) 003; 
C. F. Berger, L. Covi, S. Kraml and F. Palorini, JCAP 0810 (2008) 005. 



12 



[12] B. A. Campbell, S. Davidson, J. R. Ellis and K. A. Olive, Phys. Lett. B 256, 484 
(1991); W. Fischler, G. F. Giudice, R. G. Leigh and S. Paban, Phys. Lett. B 258, 45 
(1991); H. K. Dreiner and G. G. Ross, Nucl. Phys. B 410, 188 (1993). 

[13] W. Buchmuller, L. Covi, K. Hamaguchi, A. Ibarra and T. Yanagida, JHEP 0703, 037 
(2007). 

[14] F. Takayama and M. Yamaguchi, Phys. Lett. B 485, 388 (2000). 

[15] G. Bertone, W. Buchmuller, L. Covi and A. Ibarra, JCAP 0711 (2007) 003. 

[16] A. Ibarra and D. Tran, Phys. Rev. Lett. 100, 061301 (2008). 

[17] K. Ishiwata, S. Matsumoto and T. Moroi. larXiv:0805.1133l [hep-ph]. 

[18] A. Ibarra and D. Tran, JCAP 0807 (2008) 002. 

[19] A. W. Strong, I. V. Moskalenko and O. Reimer, Astrophys. J. 613, 956 (2004); As- 
trophys. J. 613 (2004) 962. 

[20] S. W. Barwick et al. [HEAT Collaboration], Astrophys. J. 482, L191 (1997). 

[21] B. C. Allanach, A. Dedes and H. K. Dreiner, Phys. Rev. D 69 (2004) 115002 [Erratum- 
ibid. D 72 (2005) 079902]; R. Barbier et al, Phys. Rept. 420, 1 (2005). 

[22] K. Ishiwata, T. Ito and T. Moroi. la7Xrv:0807.0975l [hep-ph]. 

[23] X. Ji, R. N. Mohapatra, S. Nussinov and Y. Zhang. larXiv: 0808 .19041 [hep-ph]. 

[24] S. Davidson and A. Ibarra, Phys. Lett. B 535 (2002) 25. 

[25] K. Hamaguchi, H. Murayama and T. Yanagida, Phys. Rev. D 65, 043512 (2002). 

[26] L. Covi, E. Roulet and F. Vissani, Phys. Lett. B 384, 169 (1996). 

[27] W. Buchmuller, P. Di Bari and M. Plumacher, Annals Phys. 315, 305 (2005). 

[28] G. F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia, Nucl. Phys. B 685, 
89 (2004). 

[29] G. Hinshaw et al. [WMAP Collaboration], larXiy^0A0732| [astro-ph]. 

[30] S. Blanchet and P. Di Bari, JCAP 0703 (2007) 018 | arXiv:hep-ph/0607330| . 

[31] M. Bolz, A. Brandenburg and W. Buchmuller, Nucl. Phys. B 606, 518 (2001) 
[Erratum-ibid. B 790, 336 (2008)]. 

[32] J. Pradler and F. D. Steffen, Phys. Rev. D 75, 023509 (2007). 

[33] V. S. Rychkov and A. Strumia, Phys. Rev. D 75, 075011 (2007). 

13 



[34] M. Endo, F. Takahashi and T. T. Yanagida, Phys. Rev. D 76 (2007) 083509. 

[35] L. Randall and R. Sundrum, Nucl. Phys. B 557, 79 (1999); G. F. Giudice, M. A. Luty, 
H. Murayama and R. Rattazzi, JHEP 9812, 027 (1998). 

[36] K. Choi, A. Falkowski, H.P Nilles and M. Olechowski, Nucl. Phys. B718 (2005) 113. 

[37] M. Endo, M. Yamaguchi and K. Yoshioka, Phys. Rev. D72 (2005) 015004; K. Choi, 
K.S. Jeong and K.i. Okumura, JHEP 0509 (2005) 039. 

[38] M. Endo, K. Hamaguchi and F. Takahashi, Phys. Rev. Lett. 96 (2006) 211301; 
S. Nakamura and M. Yamaguchi, Phys. Lett. B 638 (2006) 389. 

[39] B. C. Allanach, Comput. Phys. Commun. 143, 305 (2002). 

[40] C. Amsler et al. (Particle Data Group), Physics Letters B667, 1 (2008) 

[41] M. Frank, T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak and G. Weiglein, JHEP 
0702, 047 (2007); G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich and G. Weiglein, 
Eur. Phys. J. C 28, 133 (2003); S. Heinemeyer, W. Hollik and G. Weiglein, Eur. Phys. 
J. C 9, 343 (1999); Comput. Phys. Commun. 124, 76 (2000). 

[42] G. Degrassi, P. Gambino and P. Slavich. larXiv:0712.3265l [hep-ph]. 

[43] G. W. Bennett et al. [Muon G-2 Collaboration], Phys. Rev. D 73, 072003 (2006). 

[44] K. Hagiwara, A. D. Martin, D. Nomura and T. Teubner, Phys. Lett. B 649, 173 
(2007). 

[45] For a recent discussion and references, see M. Passera, W. J. Marciano and A. Sirlin, 
larXiv:0804. 11421 [hep-ph]; iaxXiv:0809.4"062l [hep-ph]. 

[46] J.L. Lopez, D.V. Nanopoulos and X. Wang, Phys. Rev. D 49, 366 (1994); U. Chat- 
topadhyay and P. Nath, Phys. Rev. D 53, 1648 (1996); T. Moroi, Phys. Rev. D 53, 
6565 (1996); Erratum-ibid, D 56, 4424 (1997). 

[47] W. Buchmuller, K. Hamaguchi, O. Lebedev and M. Ratz, Nucl. Phys. B 699 (2004) 
292. 



14 



